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LIBRARY ROUTINE V h - 178 



TITLE 

TYPE 

DURATION 

LENGTH 

ACCURACY 

TEMPORARY STORAGE 

ENTRY 



Fourier Analysis 

Closed 

(m + 13m) milliseconds 

52 words 

3 x 10" 12 

through 6 and m consecutive locations specified by preset 

parameter S6. 

When this routine is located at q entry is made by 



PRESET PARAMETERS 



p 


XO mF 
50 pF 


p+1 


26 qF 



where X = 5 for sine series; X = J for cosine series. 
S3 00. F If the function to be analyzed is f(x), 

00 aF < x < L, values f . = f (iL/m) > < i < m, 
of the function are stored at consecutive 
locations beginning with a. (There are m + 1 
values of the function at equal intervals 
of length L/m along the interval < x < L. 
The first value is f(0), stored at a. The 
last value is f (it), stored atm + a.) 
00 F This routine uses library routine T 5 - 157 
00 bF to compute sine or cosine values, b is 

the location of the first word of routine T 5 
00 F This routine computes values of the co- 
00 cF efficients of the sine or cosine expansion 
of a function. Scaled values of these 
coefficients are stored by the routine 
at consecutive locations beginning with c. 
For sine series the number of coefficients 
computed is the greatest integer less than 
or equal to m/2 . One more than this is 



Sk 



S5 



-2- 



computed for cosine series. 
S6 00 F This routine prepares and uses a table 
00 dF of sine or cosine values. This table 
occupies m locations beginning with d. 
NOTES ON USE This routine computes the first approximately m/2 

coefficients of the sine or cosine expansion of a function. 



oo 



f(x) = E A sin 2*£ (0<x<L) (l 
n=l n L 



or 



f(x) = B- + I B cos ~ (0 < x < L) (2 
, n L ' K 

n=l 

Values of the function f(x) at equal intervals along 
the interval from to L are stored prior to entering 
the routine. 

For sine series the coefficients are scaled by l/2 . 
For cosine series all coefficients except the first, B n , 

are scaled by l/2 . B^ is unsealed. 

u 

The accuracy of the computed values of the coefficients 
depends upon the function being analysed. The sine 
series represents the function -f(-x), and the cosine 
series the function f(-x), in the open interval (-L, 0). 
The sine series represents a function which is zero 
at the end. points, and L, and the cosine represents 
a function whose first derivative is zero at these points. 
Both series represent functions which are periodic of 
period 2L. The accuracy of the computed values of 
the coefficients depends upon the continuity of this 
function and the value of m chosen. 

More precisely, the error in the integration 
formula used (the trapezoidal rule) is 



-3- 

E 

i=l 



AA n = t (A 21m + n - Ata . ft)' (0<n<VD 



A B o " S n *2ta 

1=1 

oo 

A B_ = S (B OJ „ ^ „ + BL 4 „ J (0 < n< m/2) 



n 



1=1 



2im + n 2im ■„ n' 



Thus the error tends to "be large if the higher (relative 

to m) harmonics are not negligible. For example* the 

relative error in AV for a function whose sine coefficients 

fall off as l/n is about 0-5$ with m = 32 1 whereas 

if the sine coefficients fall off as l/n the relative 

-error in A, is about 0.0015$ for m = 32. It will also 

be noticed from the error formula that if, for instance/ 

the even harmonics of the function are zero, the error 

due to the integration formula for their even coefficients 

will also be zero. The only error for these coefficients 

in this case will be the error due to the routine 

-12 
itself, which has a maximum value of 3 x .10 

A complete Fourier series representation of a 
function, 

oo 

f (x) = B Q + Z [A n sin ^_ + B n cos — ] 
n=l 

may be obtained using this routine. The sine coefficients 
are obtained by finding the sine coefficients of the 
function 

f (x) = 1/2 [f(x) - f (-x)] (0 < x < L) 

The cosine coefficients are obtained by finding the 
cosine coefficients of the function 



f (x) = 1/2 [f (x) + f (-x)] (0 < x < L) 



The series thus obtained will represent the function 



nitx . „ __„ mtx_(-L < x < L) 



-it- 



over the interval (-L < x < L.) and will be periodic 
of period 2L. However, to obtain this series, values 
of the function must be known at an odd number (2m + l) 
of points, m itself may be either even or odd. 

This routine leaves with the address, specified 
by preset parameter S5* of the first coefficient 
computed as the right hand address of the word in 
the accumulator. 



DATE April k, 1955 



CODED B Y QA^fT^r^^p 
APPROVED BY ytipflod/s 



LOCATION 


> ~ ' ' ■ — 
ORDER 


■ — 


NOTES 


PAGE 1 





4l 3F 
4l 4f 








1 


4l 5F 
L5 49L 








2 


40 10L 
46 5F 




Set n = 1 




3 


K5 F 
42 47L 








4 


46 4F 
10 20F 








5 


42 3F 
36 8L 




Sine or cosine? 




6 


19 3F 
L4 10L 


• 






7 


40 10L 

41 5F 




For cosine values set 


n = 


8 


4l 6F 


from 5 








L5 6F 


from 15 


Prepare table 




9 


JO 6f 
66 3F 








10 


S(5) IF 

50 10L 


*>y 2, 7 






11 


26 s4 
4o ( )s6 


ty 13, 16 






12 


F5 6f 
4o 6f 








13 


F5 11L 
42 11L 








14 


L5 3F 
F0 6f 








15 


32 8L 

L5 51L , 








16 


42 11L 
L4 4f 








17 


46 37L 










L5 5F 


from 46 


Compute coefficient 





V 4 



LOCATION 
18 

19 
20 
21 
22 
23 
2k 

'25 

26 

27 
28 

29 

30 

31 
32 
33 
3^ 
35 



ORDER 

fH 20L 
Ho 2^L 
F5 50L 
kO 6F 
50 S6 
7J S3 
HO 2F 
50 2F 
75 50L 
H-0 IF 
S5 F 
kO F 
50 ( )S6 

7J ( )S3 

HO 2F 
50 2F 
L5 F 
7^ 6F 
Ln IF 
HO IF 
kl 2F 
L5 2^L 
F*J- 5F 
40 2^L 
LO 48L 
k2 2F 
LO HF 
36 3^L 
Lk kQL 
k6 2^L 
LI 6f 
40 6F 
L5 3F 
FO 2F 
36 23L 
S5 S5 



NOTES 



PAGE 2 



V k 



Integration 



from 35 

by 32 

by 18, 29 



from 31 



LOCATION 


ORDER 


NOTES PAGE 3 V 4 


36 


40 F 
50 S6 






■■ 37 


75 ( )S3 
10 IF 


by 17 




38 


40 2F 
50 2F 






39 


L5 F 
74 6F 






40 


L4 IF 








66 3F 






41 


S5 F 








40 ( )S5 


by 42, 47 




42 


F5 41L 
42 41L 




Prepare to compute next coefficient 


43 


L5 5F 

L4 49L 






44 


46 5F 
L5 4F 






45 


10 IF 
LO 5F 






46 ■ 


32 17L 
L5 35L 






47 


42 4lL 








22 ( )F 


by 3 


Exit 


48 


JO S6 
7J S3 




Test constant 


49 


S5 IF 
50 10L 




Order used for sine series only 


50 


00 F 
00 IF 




Scaling constant 


51 

s 


00 S3 
00 S6 







